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The multi-orbital Hubbard model in one dimension is studied using the numerical diagonal- 
ization method. Due to the effect of the crystal-field splitting A, the fully polarized ferromag- 
netism which is observed in the strong coupling regime becomes unstable against the partially 
polarized ferromagnetism when the Hund's rule coupling J is smaller than a certain critical 
value of order of A. In the vicinity of the partially polarized ferromagnetism, the orbital fluc- 
tuation develops due to the competition between the Hund's rule coupling and the crystal-field 
splitting. The superconducting phase with the Luttinger liquid parameter Kp > 1 is observed 
for the singlet ground state in this region. 
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Recently, the strongly correlated electron systems with 
orbital degrees of freedom have attracted much inter- 
est. Striking phenomena of such systems are the colos- 
sal magnetoresistance in manganites Lai_a:Sra;Mn03,^'' 
the triplet pairing superconductivity in the ruthenate 
Sr2Ru04^^ and the metal-insulator transition in alkali- 
doped fuUerides A^Ceo-^^ 

Another interesting class of materials is cobalt oxides 
such as Lai_a;Sra;Co03^"^-' which shows the spin state 
transition. Because the Hund's rule coupling and the 
crystal-field splitting between the t2g and eg orbitals are 
close to each other, the spin state of the cobalt ion de- 
pends on temperature, doping concentration and crystal 
structure. For example, the ground state of the Co'^+ 
(3d^) ion in LaCoOa is known to be a low-spin (LS) 
state (iigSgj ^ = 0), while, with increasing temperature, 
the spin state gradually changes into an intermediate- 
spin (IS) state (^2363' 5* = 1) and/or a high-spin (HS) 
state (4ge2, 5* = 2).5-8) With Sr doping, Lai-^Sr^CoOa 
shows a spin-glass for x < 0.18 and a ferromagnetism for 
X > 0.18,^-' where IS and/or HS states play crucial roles. 

The recent discovery of large thermoelectric power 
in Nao.5Co02"'^°'' has stimulated considerable attention 
on layered cobalt oxides such as Naa;Co02. Koshibae 
et a/.^^^ claimed that large degeneracy of electronic 
states due to a competition between the Hund's rule 
coupling and the crystal-field splitting, together with 
strong correlation plays a key role in Na2,Co02 as well 
as Lai_2:Sr2:Co03. Weak ferromagnetism has been ob- 
served in Nao.75Co02."'^^'' More recently, Takada et al}"^^ 
have discovered the superconductivity in Na2;Co02 • 
2/H2O with Tc w 5K for x « 0.35 and y « 1.3. With 
the advent of the new findings, theoretical studies of the 
interplay of the Hund's rule coupling and the crystal-field 
splitting is highly desirable. 

The orbitally degenerate Hubbard model has been 
extensively investigated to clarify the effect of orbital 
degrees of freedom in the presence of the intra-atomic 
Coulomb interaction. Many authors^^"^^-* have studied 



the ferromagnetism of this model and revealed that the 
Hund's rule coupling plays a crucial role in the ferro- 
magnetism; however the effect of the crystal-field split- 
ting was not considered there. Possible mechanisms of 
superconductivity have been proposed by several au- 
thors, ^'^"^'^^ but the relationship with the ferromagnetism 
was not discussed there. 

In the present work, we investigate the ferromagnetism 
and the superconductivity in the multi-orbital Hubbard 
model, in particular paying attention to the effect of the 
interplay of the Hund's rule coupling J and the crystal- 
field splitting A. As the strong correlation effect plays 
crucial roles in ferromagnetism, a non-perturbative and 
reliable approach is required. We employ the numeri- 
cal diagonalization method for the one-dimensional Hub- 
bard model with finite system sizes. This approach has 
already been applied for the A = case.^**' Although 
the available system size is fairly small, the results are 
in good agreement with the strong coupling analysis^^-* 
and the results from the density-matrix renormalization- 
group method. To examine the superconductivity, we 
calculate the critical exponent of the correlation func- 
tions Kp based on the Luttinger liquid theory.^'*' The 
reliability of this approach has been extensively tested 
for various one-dimensional models such as the Hubbard 
model,26) the t-J model,^^) the d-p model,2^"3°) etc. We 
can thus expect that this approach is reliable for multi- 
orbital Hubbard model as well. 

We consider the following Hamiltonian for the one- 
dimensional multi-orbital Hubbard model: 
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H (4m,<TC^+l,m,^ + h.C.) 
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- J' Yl i4,a,A,a,l^iMCi,b,l + h-C.) (1) 

where c- „ ^ stands for creation operator of a electron 
with spin a in the orbital m (= a,h) at site i and 
ni,m,<T = cJ,m,<jCi,m,o-- Here, i represents the hopping 
integral between the same orbitals and wc set t = 1 in 
this study. The interaction parameters U, U' , J and J' 
stand the intra- and inter-orbital direct Coulomb interac- 
tions, the exchange (Hund's rule) coupling and the pair- 
transfer, respectively. A stands the energy difference 
between the two atomic orbitals, i.e., the crystal- field 
splitting. For simplicity, we impose the relations, J = J' 
and U = U' + 2 J, which holds exactly in 3rf-orbitals for 
A = and is a good approximation for A ^ 0. The 
model eq.(l) is schematically represented by Fig. 1(a). 

To carry out a systematic calculation, we use the pe- 
riodic boundary condition for A'^,. = 4m + 2 and the an- 
tiperiodic boundary condition for Ng = 4to, where A^e 
is the total electron number and m is an integer. This 
choice of the boundary condition removes ac;cidc;ntal de- 
generacy so that the ground state might always be a sin- 
glet with zero momentum except for the ferromagnetic 
state at large J. The filling n is defined by n = N^/Nu, 
where Nu is the total number of unit cells (each unit cell 
contains two orbitals). We numerically diagonalize the 
Ifamiltonian eq.(l) up to 14 sites (7 unit cells) using the 
standard Lanczos algorithm. 

When the charge gap vanishes in the thermodynamic 
limit, the uniform charge susceptibility Xc is obtained 

from Xc = _Eo ( + 2 , Af„ ) + Bo ( We - 2 , ) - 2 _Eo ( A^s , ) ' '^^^'"^ 

Eo{Ne, Nu) is the ground state energy of a system with 
Nu unit cells and Ne electrons. In the Luttinger liq- 
uid theory, some relations have been established as uni- 
versal relations in one-dimensional models.^'*' In the 
model which is isotropic in spin space, the critical expo- 
nents of various types of correlation functions are deter- 
mined by a single parameter Kp. It is predicted that 
the superconducting (SC) correlation is dominant for 

Kp > 1 (the correlation function decays as ^ r ^^f" ) , 
whereas the CDW or SDW correlations are dominant for 
Kp < 1 (the correlation functions decay as ~ r~^^+^''^) 
in the Tomonaga-Luttinger regime. The critical ex- 
ponent Kp is related to the charge susceptibility Xc 
and the Drude weight D hy Kp = i(7rXc-D)^/^, with 

D = where Eo{(j)) is the total energy of the 

ground state as a fimction of a magnetic flux 7V„(/).^^'' 

In the noninteracting case {U — U' = J = 0), the 
Hamiltonian eq.(l) yields a dispersion relation €^{k) = 
— 2tcos(fc) ± -J, where k is the wave vector and e+(fc) 
(e~(fc)) represents the upper (lower) band energy. When 
the lowest energy of the upper band, e+(0), is larger than 
the Fermi energy Ekp , electrons occupy only the lower 
band with fcjr = This band structure is schematically 
represented by Fig. 1(b). Hereafter, we mainly consider 
the case with £"'"(0) > E^p- 

Figure 2 shows the value ofKp as a function of J(= U') 
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(b) 

Fig. 1. Schematic diagrams of (a) the model Hamiltonian and (b) 
the band structure in the noninteracting case. 
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Fig. 2. Kp as a function of J(= U') at A = 3.5,4,4.5and 5. The 
broken line represents a weak coupling estimation for Kp . Inset 
shows the energy difference Eo{rj>) — Eo{0} as a function of an 
external flux </>. 



for several values of A at the electron density n = 8/6. 
As J increases Kp decreases for small J, while it increases 
for large J and, then, becomes larger than unity. In the 
region Kp > 1, the SC correlation is expected to be most 
dominant compared to the CDW and SDW correlations. 
When J is larger than a certain critical value, the ground 
state changes into the partially ferromagnetic state with 
S=l from the singlet state. 

In Fig. 2 the broken line represents a weak coupling es- 
timation for Kp, where Eq is calculated within the first 
order perturbation (Hartree approximation).''^' "'^^ This 
approximation shows a good agreement with the numeri- 
cal result in the weak coupling regime, ensuring the small 
finit{^-size effect of the numerical calculation. In this ap- 
proximation, the ground state energy is obtained by the 
expectation value of the Hamiltonian eq.(l) using the 
noninteracting ground state where the lower band is ex- 
clusively occupied as shown in Fig. 1(b). Then the effect 
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Fig. 3. Phase diagram of the superconducting state with Kp > 1 
(shadowed region) and the ferromagnetic state with 5 ^ on 
the U' — J parameter plane with contour map of Kp. 



of the upper band is omitted and the result is indepen- 
dent of A as shown in Fig. 2. This approximation breaks 
down when the effect of the upper band becomes crucial. 
In this regime, Kp rapidly increases with increasing J, 
and finally becomes larger than unity showing the super- 
conducting state. 

The critical values of J with Kp — 1 are Jc ~ 0.5, 
0.9, 1.3 and 1.6 for A ~ 3.5, 4, 4.5 and 5, respec- 
tively. We find that the phenomenological equation 
Jc = e^(0) — Ekp can approximately lead to the above 
results. Here, €+(0) — Ekp corresponds to the lowest 
energy of the single-particle excitation from the lower 
band to the upper band. This equation is found to be 
a good approximation for various filling n. When J ex- 
ceeds e+(0) — Ekp, the inter-band excitation develops 
and the occupation number of the upper band increases, 
which results in the large orbital fluctuation accompa- 
nied by the fluctuation between the low-spin and the 
high-spin states. The mechanism of the superconductiv- 
ity is related to this orbital fluctuation. 

We note that the superconductivity is also observed 
for the J' = (J 7^ 0) case in contrast to the previous 
study^^'^^' where the pair-transfer J' is crucial for the 
superconductivity. The importance of the upper band 
has also been pointed out in different types of multi-band 
models such as the d-p modeP*"'^"''^^^ and the Hubbard 
ladder model.^^'^^^ 

To confirm the superconducting state, we calculate the 
lowest energy of the singlet state i?o(0) as a function of 
an external flux </>. As shown in the inset of Fig. 2, the 
anomalous flux quantization occurs at J = 1.1, where 
Kp is about 1.4. When J — 0.8 and 0.4 , Kp is less than 
unity and the anomalous flux quantization is not found. 

In Fig. 3, we show the phase diagram of the supercon- 
ducting state with Kp > 1 together with the ferromag- 
netic state with S ^ on the U' vs. J parameter plane, 
where the value of Kp is given by the contour map. The 
superconducting phase appears near the partially polar- 
ized ferromagnetic region. It extends from the attractive 
region with J < and [/' < to the realistic param- 




Fig. 4. Global phase diagram on the U' — J plane. 



eter region for 3d transition-metals with U' > J > 0. 
We have confirmed that the superconducting region in- 
creases as A decreases as shown in Fig. 2. 

Fig. 4 shows the global phase diagram on the U' — J 
plane. When U' ^ J ^ A, the fully polarized ferro- 
magnetism with iS* — S'max appears. It accompanies the 
partially polarized ferromagnetism with < 5 < Smax 
for J ^ A. In the A = case, the ferromagnetism of 
the degenerate Hubbard model has been previously stud- 
ied. In one dimension, some rigorous results are 
shown in the strong coupling limit [/ — > cxd: the ground 
state is fully polarized ferromagnetism for < n < 2 
except for n = 1 when U' and J = J' are positive and 
finite. Numerical result suggests that the ferromag- 
netism is stable also for n = 1 in the strong coupling 
region. In infinite dimensions, the dynamical mean- 
field theory shows the existence of the ferromagnetism in 
the same parameter region observed in one dimension. 
The ferromagnetism for n ^ 1 is found to be metal- 
lic and mainly caused by the double-exchange mecha- 
nism,^^) which is also the present case with A ^ 0. The 
existence of the partially polarized ferromagnetism for 
A > has been reported in a different type of the two- 
band Hubbard model. '^^^ The ferromagnetic phase for 
[/' is complicated as shown in Fig. 4. The origin of 
this phase is not clear at this stage. 

In Fig. 5, we show the doping dependence of the critical 
values of J for the superconductivity and the ferromag- 
netism at A = 4 with J — U', where we use n=8/6, 10/7, 
6/4, 8/5 and 12/7 systems. The critical values are de- 
fined above which Kp > 1 for the superconductivity and 
S > for the ferromagnetism, respectively. Although the 
finite size effect is considerably large, the phase boundary 
for the superconductivity can be approximately given by 
the phenomenological equation, Jc — e^(0) — Ekp, men- 
tioned before. 

For n = %ll and n = 6/5, we cannot use the expres- 
sion Kp = ^{nXcDy^^ , because the system has a charge 
gap at Ne — Nu and, then, Xc is not obtained from the 
finite size systems. Instead we use another expression, 
Kp — with the charge velocity Vc- Although there 
is a larger finite size effect, we find Kp < 1 until the 
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Fig. 5. Phase diagram on the n — J plane. The broken hne rep- 
resents the phenomenological equation Jc = £"''(0) — E^-p (see 
text) as a function of n. 



ferromagnetic transition occurs. Furthermore, the result 
of the flux quantization also shows no sign of the super- 
conductivity in this case. These results suggest that the 
superconducting phase may disappear for n < 6/5. 

In conclusion, we have obtained the phase diagram of 
the one-dimensional Hubbard model with two-fold or- 
bital degeneracy. The fully polarized ferromagnetism 
has been found in the strong coupling regime with U' Z 
J ^ A. For 1 < n < 2, the ferromagnetism is metal- 
lic and mainly caused by the double-exchange mecha- 
nism. The crystal-field splitting destroys the fully 
polarized ferromagnetism resulting in the partially po- 
larized one for J < A. In the vicinity of the partially 
polarized ferromagnetism, we have found the supercon- 
ducting phase, when J exceeds the lowest energy of the 
inter-band excitation, which extends to the realistic pa- 
rameter region for 3d transition-metals with U' > J. 

Sakamoto et al}^'' claimed that the metallic ferro- 
magnetism appears in a similar parameter region in any 
dimension by comparing the results from one dimension 
with those from infinite dimensions. It is natural to think 
that this ferromagnetism will appear in two and three 
dimensions even in the presence of A. Then, we expect 
that a partially polarized (weak) ferromagnetism appears 
in real materials, in which the Hund's rule coupling and 
the crystal-field splitting compete to each other, such as 
cobalt oxides. In fact a weak ferromagnetism has been 
observed in the layered Nao.75Co02^^-' as well as the per- 
ovskite Ri_:,A^Co03.^^^ 

Finally we make some comments on the supercon- 
ductivity in Naa;Co02 • yll20}^^ Although the several 
authors'^^"'^*' have discussed the system using the non- 
degenerate t-J model, the orbital degeneracy of 3d elec- 
trons is considered to play a crucial role in Naa;Co02 as 
well as Lai_2:Sra;Co03.^^'' The competition between the 
Hund's rule coupling and the crystal-field splitting causes 
the large orbital fluctuation, accompanied by the fluctu- 
ation between the low-spin and the high-spin states at 
each Co ion, which mediates the superconductivity. As 
the orbital fluctuation has a local character, the mech- 



anism for the superconductivity could be common in 
all dimensions. We hence expect that the multi-orbital 
mechanism is responsible for the superconductivity in 
layered Na2:Co02 • 2/H2O. Exploration of the supercon- 
ductivity in the vicinity of the weak ferromagnetism in 
the perovskite Ri-jjAjjCoOs^^^ may be promising. 
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cation, Science, Sports and Culture. 
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